RADIAL DUNKL PROCESSES : EXISTENCE AND 
UNIQUENESS, HITTING TIME, BETA PROCESSES AND 

RANDOM MATRICES 
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Abstract. We begin with the study of some properties of the radial Dunkl 
process associated to a reduced root system R. It is shown that this diffusion 
is the unique strong solution for alH > of a SDE with singular drift. Then, 
we study Tq, the first hitting time of the positive Weyl chamber : we prove, via 
stochastic calculus, a result already obtained by Chybiryakov on the finiteness 
of To . The second and new part deals with the law of Tq for which we compute 
the tail distribution, as well as some insight via stochastic calculus on how root 
systems are connected with eigenvalues of standard matrix-valued processes. 
This gives rise to the so-called /3-processes. The ultraspherical /3-Jacobi case 
still involves a reduced root system while the general case is closely connected 
to a non reduced one. This process lives in a convex bounded domain known 
as principal Weyl alcove and the strong uniqueness result remains valid. The 
last part deals with the first hitting time of the alcove's boundary and the semi 
group density which enables us to answer some open questions. 

1. Preliminaries 

We begin by pointing out some facts on root systems and radial Dunkl processes. 
We refer to [38] for the Dunkl theory, to both [7] and [27] for a background on 
root systems and [H], [21] for facts on radial Dunkl processes. Let (V, <, >) be a 
finite real Euclidean space of dimension m. A reduced root system i? is a finite set 
of non zero vectors spanning V such that : 

1 RnRa = {a, -a} for all a e R. 

2 aa{R) = R 

where cTq is the reflection with respect to the hyperplane Ha orthogonal to a : 

o'a[x) = X — 2 a, X E V 

< a,a > 

A simple system A is a basis of V which induces a total ordering in R. A root a is 
positive if it is a positive linear combination of elements of A. The set of positive 
roots is called a positive subsystem and is denoted by Note that the choice of 
A is not unique and that R+ is uniquely determined by A. The reflection group 
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W is the one generated by all the reflections cTq, for a & R. Recall that W is finite 
and the only reflections are of the form for a & R. Given a root system R with 
associated positive subsytem let C be the positive Weyl chamber defined by : 

C := {x eV < a,x » OVa G = {x e V < a,x » OVa G A} 

and C its closure. One of the most important properties is that the convex cone 
C is a fundamental domain, that is each A G V is conjugate to one and only one 
/i G _ 

The radial Dunkl process is defined as the C-valued continuous paths Markov 
process whose generator is given by : 

^u{x) = -An(x) + ^^^^ 

with boundary conditions Vm(x) ■ a = for all x G Ha, a G R+, k{a) > is the 
multiplicity function (invariant under the action of W), and u G C^{C). When 
k{a) = 1 for all a G -R, we recover the BM constrained to stay in C, studied by 
Grabiner (1^2^). The semi-group density of X is given by : 

(1) pU^, y) = (X, y)ll<a,y >2^(") 



for x,y eC, where 7 = EaeJ?+ ^( 



a 



X wy 



where Dk denotes the Dunkl kernel and Ck is given by the Macdonald-Mehta inte- 
gral (PI). Indeed, as Dk{0,y) = 1 ([38j), one gets 

t7+W2c^ = f ^~\y?l2t rr < >2^H dy= \ e-\y\"l^' n \<a,y> ^'^("^d; 

since = U^i^wwC . D^{x,y) is known as the generalized Bess el function (up 
to the constant \W\). This process is obtained by projecting the Dunkl process 
valued in (which has right- continuous and left-limits paths, see [2Tj) on C. The 
latter was already introduced by Rosier ([38], [39]) and then studied by Gallardo 
and Yor ([2I],[22]) and Chybiryakov ([IT]). 

To illustrate all these facts and motivate the reader as well, we will provide 
some well known examples. We start with the rank one case (m = 1) for which 
R = Bi = {±1}. Hence k{a) := ko > and X is a Bessel process ([3?]) of index 
u = ko — 1/2. When /cq > 0, it is the unique strong solution of : 

k 

dXt = dBt + -^dt, t > 0, Xo = X > 0. 



xr' 
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where i? is a standard BM. Another well known multivariate example is described 
by the A^-i-type root system defined as : 

Am~i = {±{ei - ej) 1 < i < j < m}, 

with positive and simple systems given by : 

R+ = {ci — Cj , 1 < i < j < m} A = {cj — Cj+i, I < i < m} 

where (ej)i<j<m is the standard basis of M™. V is the hyperplane of R"^ consisting 
of vectors that coordinates sum to zero. Without loss of generality, one can take 
M."^ instead of V so that C = {x E M™, xi > ■ ■ ■ > Xm}- Besides, there is only one 
orbit and k{a) := ki > 0. Thus, the corresponding radial Dunkl process satisfies : 

(2) dXf = duf + A;i V — — . 1< i < m, t<T 

with Xq > ■ ■ ■ > X™, where (z/*)j are independent Brownian motions and r is 
the first collision time . For strictly positive ki, this process was deeply studied 
by Cepa and Lepingle ([8], [9], [10]): it behaves as m-interacting particles on the 
real line with electrostatic repulsions proportional to the inverse of the distance 
separating them. Moreover, when ki = 1,1/2 respectively, this process evolves like 
eigenvalues process of Hermitian (Dyson model) and symmetric Brownian motions 
([20], [23])- It was shown in [B] that this SDE has a unique strong solution for all 
t > and ki > 0. When reading the proof in [S], one hopes to extend this result 
for any root system since materials used there are not typical for the A^-i-type. 
This was the original motivation of this work. Our first result claims that 

dXt = dBt - V^{Xt)dt, XoeC 

where $(x) = — J2a&R+ ^(^) lii(< a, x >), k > , has a unique strong solution for 
alH > 0. At the same time and independently, Chybiryakov and Schapira provide 
two other proofs: both authors used well posed martingale problems associated 
respectively with the M'^-valued Dunkl and the radial Heckman-Opdam processes 
as well as geometric arguments ([H], ^^)■ The curious reader will wonder what 
happens if k{a) = for some a? The answer is that the same result holds up to the 
first hitting time of dC, say Tq ([11] p. 37). Next, we are mainly interested in the 
tail distribution of Tq. Before proceeding, we reprove via stochastic calculus that 
To < oo if k{a) < 1/2 for at least one a G R+ (see [TT] for the original proof using 
local martingales). More precisely, for such an a, we prove that almost surely, 
< a,Xt >< Yt for all t>0, where F is a Bessel process of dimension strictly less 
than 2. At this level, other proofs exist for the above results. To our knowledge, 
the contents of the remainder of the paper are new. In [TT], the author derived 
absolute-continuity relations which allow us to write the tail distribution of Tq 
when starting from x & C. A IV-invariant analytic x-dependent integral, which 
value at is given by a Selberg integral, is involved. As far as we know, though 
D^{x,y) arises as hypergeometric functions for particular root systems (see the 
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end of [2]), forward computations are sophisticated and hard. More precisely, we 
think that it is possible to use the integral formula given in Corollary 2 in [28] 
with the integration range < < ■ ■ ■ < < 1, known as the Macdonald's 
conjecture, then perform limit and sums operations. The matrix cases for which 
the Jack parameter equals to = 1,2 are more handable with the use of properties 
of zonal polynomials and Schur functions. However, we think that the approach 
adopted here is more elegant since on one hand, it disgards the special values of 
the multiplicity function and on the other hand, does not need long hard formulas. 
It only relies on some properties picked from Dunkl theory. More precisely, it will 
be shown that the x-dependent integral is an eigenfunction of some operator which 
involves the generator ^ and the so-called Euler operator Ei. For some particular 
root systems, this eigenfunction is identified with some hypergeometric series. A 
surprising fact is that the eigenoperator can be expressed in terms of a Schrodinger 
operator and its minimal eigenvalue Emin (minimal energy) (see [3H] page 18): 



1=1 

Moreover, {Xt)t>o specializes for some values of k to eigenvalues processes of self- 
adjoint matrix processes such as symmetric and Hermitian Brownian motions, 
Wishart and Laguerre and matrix Jacobi processes. In those cases, computations 
can be performed using the action of orthogonal and unitary groups. Indeed, Jack 
polynomials fit zonal polynomials and Schur functions when the Jack parameter 
equals to 1, 2 respectively (see [34]). The two first ones are identified as A^-i-type 
radial Dunkl processes while Wishart and Laguerre processes are related to the 
Bm root system. The latter goes beyond the radial Dunkl setting: the reduced 
root system Cm in a particular case (ultraspheric) is involved and more generally, 
the non reduced system BCm- This connection was deeply investigated in [3] 
while identifying special functions associated with root systems with multivariate 
hypergeometric series. Among them appear multivariate Gauss hypergeometric 
series and Jacobi polynomials ([S2]) and these are eigenfunctions of the /3- Jacobi 
generator. The state space is the so-called principal Weyl alcove which is now a 
bounded convex domain and fundamental for the action of the affine Weyl group. 
Hence, the process evolves like particles in an interval. Then, we extend the 
strong uniqueness Theorem to the Jacobi context. In the remaining part, we 
derive some properties: we briefly visit the Brownian motion in the principal Weyl 
alcove which corresponds to multiplicities all equal to 1. Then, an analogous result 
on the finiteness of the first hitting time of alcoves walls is obtained using similar 
computations as those for Tq. Finally, we derive the semi group density and discuss 
some open questions left in [TB] . 
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2. Radial Dunkl Process : Existence and Uniqueness of a strong 

solution 

Theorem 1. Let R be a reduced root system. Let: 

<l>(x) = - ^ A;(a)ln(< a,x >) := ^ k{a)9{< a,x>), x eC 

where k{a) > for all a G R+. Then the SDE 

(3) dXt = dBt - V^Xt)dt, X^ e C 

where X is an adapted continuous process valued in C and B is a Brownian motion 
in M™, has a unique strong solution. 

Proof. From Theorem 2. 2 in [9j, we deduce that : 

(4) dXt = dBt - V^{Xt)dt + n{Xt)dLt, Xq e C 

where n{x) is a (unitary) inward normal vector to C at x , L is the boundary 
process defined by: 

dLt = l{Xt<=dC}dLt, 
has a unique strong solution for alH > 0. Moreover : 



(5) E 



(6) E 



T 

T 



dt 



\V<t>{Xt)\dt 











< oo 



for all T > 0. Thus, it remains to prove that the boundary process vanishes. To 
proceed, we need two Lemmas. 

Remark. Both Lemmas below discard the reducedness of R. In fact, this assump- 
tion figures in the definition of the Dunkl process and originates from analytic 
purposes like the commutativity of Dunkl operators ([18]). 

Lemma 1. Set dGt '■= n{Xt)dLt. Then, Va G R+, 

l{<Xt,a>=o} < dGt, a >= 

Proof: The proof is roughly a generalization of the one in [S] for R = A^-i ■ In 
order to convince the reader, we provide an outline. Using the occupation density 
formula, we may write (< a,X >> 0) : 

/•oo /»i 

/ L1{< a,X >)9' {a)da =< a,a > / 9'{<a,Xs>)ds 
Jo Jo 

where Lt{< a,X >) is the local time of the real continuous semimartingale < 
a,X >. On the other hand, the following inequaliy holds (instead of (2.5) in [S]) 



for all a G C: 

< V$(a;), X — a > = k{a)6 {< a,x >) < a,x — a > 

aeR+ 

> y, ^(o^) i^aO {< a, X >) — Ca < a, X — a > —da] 

> min {bak{a)) 6* (< a, a; >) — |x — a\ fc(Q;)cQ,|Q;| — 

:= A ^ 9'{< a,x>) - B\x - a\ - C 

by Cauchy-Schwarz inequality, where in (1), we used eq. (2.1) in [8] : let g be 
a convex C^-function on an open convex set D C M™, then Va G -D, there exist 
b,c,d > such that for all x G -D : 

< 'Vg{x),x — a > > b\Vg{x)\ — c\x — a\ — d 

Note also that A > since bak{a) > for all a G -R+. Then, the continuity of X, 
and the inequality above yield : 



which implies that 



/ 9'{< a,Xs >)ds < oo 
Jo 

POO 

/ L^{< a, X >)9' {a)da < oo 
Jo 



Thus, L'^{< a, X >) = since the function a 9' (a) is not integrable at 0. The 
next step consists in using Tanaka formula to compute dZt :=< a,Xf > — (< 
a,Xt >)+ for a G A : 

dZt = l{<a,Xi>=o} < a,dBt > -l{<a,Xt>=Q} < a,V^{Xt) > dt+l{^a,Xt>=o} < Oi,dGt > 

It is obvious that the second term vanishes. The first vanishes too since it is a 
continuous local martingale with null bracket (occupation density formula). As 
Xt G C, then dZt = a.s. which gives the result. ■ 

Lemma 2. Let x G dC. Then < n{x),a >^ for some a G A such that 

< x,a >= 0. 

Proof : Let us assume that < n{x), a >= for all a G A such that < x,a >= 
0. Then, our assumption implies that < x,a » for all a G A such that 

< n{x), a >7^ 0. If < n{x), a >< for these simple roots, then x — n{x) G C. By 
the virtue of the definition of the inward normal n{x) to C at x, i. e, 

(7) < X - a,n{x) ><0, 'iaeC, 
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it follows that n{x) is the null vector which is not possible. Otherwise, choosing : 

< x,a > 



< e < min 



a/<n(xO,a>>0 < n{x) , a > 

we claim that a := x — en{x) G dC. Arguing as before, we are done. ■ 
Now we proceed to end the proof of the Theorem. Note first that dC = Ua^AHa 
so that : 

'^{Xt<^dC}dLt < ^ l{^Xt,a>=Q}dLt. 

aeA 

If Xt G Ha for one and only one a G A. Then, n{Xt) = a/\\a\\ and Lemma [H 
gives 

l{<Xt,Q>=o} < dGt,a >= l{<Xt,o>=o}||«IM-^t = 

Hence, the boundary process vanishes. More generally, we can use the inequality 
above and write 

< Lt / l{<Xs,a>=0}'i-{<n{Xs),a>^0}dLs 



/ l{<n(X«),a>^0}— TYY^— l{<X„a>=0} < dG„a >= 



by Lemma [H 



Remark. When m = 1, {Xt)t>o is a Bessel process of dimension S = 2fco + 1 and 
fco > (5 > 1. It is well known that the local time vanishes (see Ch. XI in [37j ) 
which fits our result. 



3. FiNITENESS OF THE FIRST HITTING TIME OF THE WEYL CHAMBER 

Let To := inf{t > 0,Xt E dC} be the first hitting time of the Weyl chamber. 
It was shown in [TT] (see p. 30) that Tq = oo almost surely if k{a) > 1/2 for all 
a G R+. In [9], where R = Am-i and Tq = inf{t > 0, XI = Xj for some (i, j)}, 
authors showed that Tq < oo a.s. if and only if < /ci < 1/2. More generally, the 
following holds (see [llj p. 75 for the original proof) : 

Proposition 1. Let E A and Tq,,, := inf{t > 0, < aQ,Xt >= 0} such that 
Tq = infcQgA^ao- //" < k{ao) < 1/2, then (< ao,Xt >)t>o hits almost surely 0. 
In particular, Tq < T^y < oo a. s. 

Proof : assume k{a) > for all a G i? and let uq G A. Our scheme is roughly 
the same as that used in [9], thus we shall show that the process < aQ,X > is 
almost surely less than or equal to a Bessel process with dimension 2/c(ao) + 1. 
The result follows from the fact that 2A;(ao) + 1 < 2 when k{a) < 1/2. For this, 

7 



we use the SDE ©. For all t > 0, 

7 I I I I 7 7 / \ < CtO > 7 

a < ao, Xt > = I |a;o| |a7i + ^i'^) v '^^ 

a€R+ <Oi,Xt> 

= \\ao\\d^t + ko — Ll^^^^lL A;(q;) '^^'^^ 

<ao,At> ^ <a,At> 

ae/<+\ao 

where /cq is the value of k{ao) corresponding to the conjugacy class of ao- Set 

R = U^j^.R^ 

where R^ , I < j < p denote the conjugacy classes of R under the VT- action, then 

R+ = U^=iR+ 

so that: 

I I I 12 P 

7 1 1 1 1 7 7 K^o 7 7 < a,ao > 

a < ao, Xt >= WaoWd'jt + kf) — — at + } kj > — — at 

<ao,Xt> ^ ' ^ <a,Xt> 

For a conjugacy class R^ and a G R\ if < a,ao >= a{a) > then, it is easy to 
check that < ao{a),ao >= —a{a) where ctq is the reflection with respect to the 
orthogonal hyperplane defined by : 

<x,ao> 

aQ[x) = X — 2 ao 

< ao, ao > 

Note that <Jo{a) belongs to the same conjugacy class of a and that o"o(a) G R+ for 
a G R+ \ ao- Indeed, o"o(-R+ \ ao) = R+ \ ao for all ao G A (see Proposition 1. 4 
in [27]). Hence, 

7 M 117 7 ll'^olP 7 7 a(a) < a — (Jo(a) , Xt > , 
a < aQ,Xt >= WaoWa'jt+Ko — — at— > A;,- > — , , — at 

a{o)>0 

Furthermore, 

<a,ao> ^ N ^ ^ ^ / x<«o,^t> 

a — (Tolctj = 2 ao ^ < a — cro(aj, Aj >= za(aj — j-j — -tz — 

< ao, ao > I |tto| 1 

Consequently, one gets : 

1 1*^01 P 

d < ao, Xt >= \ \ao\\d'yt + 7? — dt + Ft dt 

< Oio,Xt > 

where Ff < on {Ta^ = oo}. Using the comparison Theorem in [29j (Proposition 
2. 18. p. 293 and Exercice 2. 19. p. 294), one claims that < ao,Xi >< ^^f^^ipi 
for all t > on {T^^ = oo}, where is a Bessel process defined on the same 
probability space with respect to the same Brownian motion, of dimension 2^0 + 1 



and starting at Yq = x >< ao,XQ » 
process of dimension < 2 hits a. s. ([37J 



0. This is not possible since a Bessel 
Chap. XI) ■ 



Remark. If we remove the assumption k{a) > for all a E R, then the SDE ([3]) 
can be solved up to time To when starting from x G C (see [II]). 

4. The law of Tq 

Here, we focus on the tail distribution of Tq deduced from absolute continuity 
relations derived in ([H]). Recall that (see [H]) the index of X is defined by 
I (a) := k{a) — 1/2. The last result asserts that if —1/2 < I (a) < for some 
a G A, then Tq < oo almost surely. Besides, if I (a) > for all a G A then 
Tq = oo almost surely. Taking into account these statements, two major parts 
are considered: l{a) > for all a G i? so that the process with index — / hits 
almost surely, and l{a) < for at least one a. The tail distribution involves a 
l^-invariant x-dependent integral. Our line of thinking relies on showing that it is 
an eigenfunction of an appropriate differential operator. Then, using uniqueness 
results for some differential equations, the tail distribution is written in Am-i and 
Bm cases by means of multivariate hypergeometric functions. In the last case, 
we recover known results from matrix theory for Wishart and Laguerre processes. 
However, we find it better to postpone this in the next section where links with 
eigenvalues of matrix-valued processes are detailed. 

4.1. A first formula. Let us denote by the law of (Xf)j>o starting from x G C. 
Let El be the corresponding expectation. Recall that ([H], Proposition 2.15.c), if 

l{a) > OVa G R+, then: 

-21(a)' 

< a,Xt> ' 



n 



< a,x > 



n 
n 



< a, X > 



21(a) 



]x\^/2t 



< a,x > 



21(a) 



Ckt-^^^/^ Jc 

g-|x|2/2t 



e-l^l^/^*/^r(^, ^) n <«'^>^^ 



n 



Ckf^ JC 

/3 I I / 2t / rt^ 



e-l«P/^D«'(^.y) II <a,y>dy 



where 7 = J2a&R+ ^(") and 7' = 7 - 1^+1/2- 

Though is given by hypergeometric functions in the special cases Am-i and 
Bm (see the end of [2]), the Jack polynomials defining them prevent us from making 
computations. However, this may be possible when these polynomials fit, for some 
values of k, Zonal polynomials and Schur functions (see [31] for definitions). Our 
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main result does not make these restrictions and uses some properties of the Dunkl 
kernel Z)^: 



Theorem 2. Let Ti be the i-th difference Dunkl operator and = YlltLi Tf the 
Dunkl Laplacian f|38] ). Define : 



i=l 



where Ef := Y^iLi^i^^f ^■^ the Euler operator and the superscript indicates the 
derivative action. Then 



Proof: Recall that if / is ly-invariant then T^f = dff and that T^Dk{x,y) 
yiDk{x,y) (see |38]). Then, on one hand : 



AlD^{x,y) = J2 Y.^wy){r^D,{x,wy) = J] J2^wyrM^,wy) 
wew i=i wew i=i 

m 

= ^yiYl Dk{x,wy) := p2{y)D^{x,y) 



i=l weW 



On the other hand 



E^D^{x,y) = J2 = J] J2^x,){wy)Mx,wy) 

weW i=l w€W i=l 

= ^ < x,wy> Dk{x, wy) = ^ < w~^x, y > Dk{x, wy) 
= ElDf{x,y) 



where the last equality follows from Dk{x, wy) = D^iw ^x, y) since D^iuix, wy) = 
Dk{x, y) for all w G W. The result follows from an easy computation. 

Corollary 1. g is an eigenfunction of — corresponding to the eigenvalue m + 

\R+\- 
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Proof : Theorem 1 and an integration by parts give : 



Ey\e-^y\y'D^{x,y)] J] <a,y>dy 



-J2 y U <a,y>df[e-\y^'/'D]r{x,y) 



dy 



Vt W <a,y> 



lit n 

^=l 

n lib 

•'^ a£R+ 1=1 



dy 



E 



oiiyi 



<a,y> 



dy 



and the proof ends by summing over i. ■ 

• The Am-i case : as mentioned in the introduction, the A^-i-type root 
system is characterized by : 

R = {±(ei — Cj), I < i < j < m} i?+ = {cj — ej, 1 < i < j < m} 

A = {ei- Ci+i, 1 < i < m} C = {y e W"^ , yi > y2 > ■ ■ ■ > ym} 

W = Sm is the permutations group and there is one conjugacy class so 
that A; = A;i>0^7 = kim{m — l)/2. Moreover, the generahzed Bessel 
function I is given by ([2] p. 212-214, [T2]): 



^ D^{x,y) = oFl,'^'^\x, 



W 



p=0 T 



where r = (ti, . . . , r^) is a partition of weight |r| = p and length m, J^^^^^ 
is the Jack polynomial of Jack parameter l/fcjl, (see [2j, ^34j)- Hence, 
letting V to be the Vandermonde function : 



p-'(To > t) = l^(a;)2^°-^^^ 



Besides, writes on ly-invariant functions 



e~\y\y\Fr''\^^,yMy)dy 



m ^ m 



i=l 



. . . Xi Xj 



i=l 



"^Authors use the change of variable x i— > \/2x, y i— > \/2y to fit the hypergeometric function 
obtained when deriving the generating function for Hermite polynomials. This in turn will modify 
the eigenoperator by a multiplying constant (see p. 183). 

''with the same notations in ki = 2 /a. This can be seen either from the eigenoperator 
below or from the orthogonality weight function involved in the semi group density. 
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Finally, since g is VT- invariant, then 

{Dl-El)g{x) = 



m + 1 , , 
m — - — g[x) 



9(0) 



\y\''^V{y)dy 



c 



\V{y)\dy 



Let us recall that the Gauss hypergeometric function 



P=0 T 



is the unique symmetric eigenfunction that equals to 1 at of (see [3] p. 
585) 

J2 z^{l-z,)df'+2h J2 ~ [c-h{m-l)-{e + b + l-h{m- 1)) z,] 



i=l 



. , . Zi 



i=l 



associated to the eigenvalue meb. Letting z = (1/2)(1 — x/^/b), e = (m + 
l)/2 and 

, / -, N 1 r , , , / -, M b ki , ^ , m + 3 
c = ki{m -l) + -[e + 6 + l- ki{m - 1)] = - + ^ (m - 1) + 

the resulting function is an eigenfunction of 



(1 - xf/b) 



1=1 l^J 1=1 

and Dq — Ef is the limiting operator as b tends to infinity. Hence, 
Proposition 2. For ki > 1/2, 



dt-Y.(b + ^-ki{m-l))^d: 



g{x) = giO)C{m, ki) lim 
where 



(lAi) 



X 



h ki , , m + 3 1 , 

m + 1, 6, - H (m - 1) H , - 1 ^ 

''2 2^ ^ 2 '2V ^/b 



C{m, k^)-' = lim sF^'/'^^ ( m + l,b,l + ^{m ~ 1) + 



m + 3 1 
2 '2 



• The i?m case : This root system is defined by 

R = {icj, icj ± Gj, 1 < i < j < m} i?+ = {cj, I < i < m, Ci ± Cj, 1 < i < j < m} 
A = {a- Ci+i, 1 < i < m, e„} C = {y G M™, yi > ^2 > ■ • ■ > Z/m > 0} 

The Weyl group is generated by transpositions and " change sign" reflec- 
tions (xj —Xi) and there are two conjugacy classes so that k = (fco, ^i) =^ 
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7 = mkQ + m{m — l)ki. The generalized Bessel functioE0 is given by ([2] 
p. 214) : 

— (x, y) = oFi'/''\ko + (m - l)k, + J, ^, ^) 



\W\ " ^ ^ ' ' ' 2' 2t' 2t 

where 

and (c)t- := Hi^ilc — — l))Ti is the generahzed Pochammer symbol (see 
[2]). Then, one has : 

r 1 2 2 ™ 

g{x) = \W\ / e-l^l'/Vi^^/'=^^(A:o + (m-l)fci + -,y,|-)n(l/.)ny')rfy 



The eigenoperator writes on Pi^-invariant functions: 



m m- 

= Y.9'^'' + 2koj2-d:+2k,j2 

■ 1 ■ 1 ■ I ■ 

%=\ x=\ l^J 



1 1 

+ 



X, i~f . rf . I />"» , 

m 



9f - 



1 /" 

-/^g{x) = mim+l)g{x), g{0) = e'^y^' llMV{y')\dy. 

A change of variable = y/2y'i shows that u{y) := g{^/2y) satisfies 
-jMy) = m^^!^^±^u{y), giO)=uiO) 

1=1 i^j ^' ^ ^ i=l 

which implies that : 



where 



u{y) = u(0),f',"'"'C.^, h + (m - l)h + \,y) 



p=0 T ^ ^'^ 



This can be seen from the differential equation ([8]) and using ([2j): 



■^there is an erroneous sign in one of the arguments in p]. Moreover, to recover this expression 
in the Bm case from that given in [2J, one should make substitutions a = fco~l/2, ki = 1/a, q = 
1 + (to — l)fci. We point to the reader that this is different from the one used in 12J p. 121. 
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Finally 



Hence, the tail distribution is given by : 
Proposition 3. Forko,ki > 1/2, 



i=l 



V — 

\2t 



Remark. 1 / Adopting the notations used in [2] , one has : 

-J'k = Dl + {a + l)El~El {R = B^,y = x'), 

Besides, Theorem |2] was derived there differently for both Am-i and 
cases when proving a generating function Theorem for generalized Hermite 
and Laguerre polynomials (page 183 and 192, see also [12j). 

4.2. A second formula. In [TT] (see Proposition 2.15.b), the author derived an- 
other absolute-continuity relation from which we deduce that if l{a) < for at 
least one a G -R+, then 



Pi(To>t) = E° 



e: 




-ds 



< a, 7 > /(a)Z(7) 
< a,Xs >< 7,Xs>' 

* < a, 7 > /(a, 7) 
< a,Xs >< 7,^s >' 



-ds 



where the last equality follows from part (c) of the same Proposition, /(a, 7) 
Z(a)Z(7) — r{a)r{'y) and 



l{a) if l{a) > 
-l{a) if /(a) < 



Then /(a, 7) = if l{a)l{'y) > and /(a, 7) = — 2r(a;)r(7) else. As a result, 

/ 



Pi(To>t) = i?: 



n 

i(a)<0 



< a, X > 



2r{a) 



exp 



E 



V(a)/{7)<0 



< a, 7 > r(a)r(7) 
< a,A:, >< 7,^s > 



(is 



Next, note that the exponential functional equals 1 for both root systems Am-i 
and Bm- For the first, it is obvious since R consists of one orbit so that {a, 7 G 
i?+,/(a)Z(7) < 0} is empty. This gives the same expression already considered in 
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the previous subsection. For the second, writing = {cj, 1 < i < m} U {e^ ± 
^kA ^ j < k < m} so that < Cj, Cj ± >= 6ij ± gives 



i=l i<k 

rn 



1 1 

+ 



-1 



i=l i<k 



EE 



=1 k<i 



= 



where 5* stands for the sum between brackets. The reader can also check that 
this holds for Cm and D^, root systems (see the end of the paper for definitions). 
However, we restrict ourselves to the B^, -case since, for particular values of the 
multiplicity function, we will recover a known result from matrix theory (see next 
section). Let us investigate the case < 1/2, ki > 1/2 for which Iq < 0, h > 0. 
One writes : 

P m 

g{x)= / e-\y\'/'DY{x,y)ll{y,)V''^{y')dy 

The machinery used before still applies and gives : 

-/kg = 2m[l + ki{m- l)]g 

Thus 

Proposition 4. In the B^, case and for ko < 1/2, ki > 1/2, one has: 



/ 2\ko-l/2 

Pi{T,>t)^CuY[{f^ 



In the remaining case k^ > 1/2, ki < 1/2, the tail distribution writes : 

» m 

g{x)^ / e-\y\'/'D^{x,y)Yl{y,r>^''V{y')dy 

Thus 

-c/kg = m[2kQ + m]g 

so that 

Proposition 5. For ko > 1/2, ki < 1/2, 

Yt) ^"'^''^'*i^i'^'^Hfco + f , ^0 + (m - l)k, + |) 
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5. /3-PROCESSES AND RANDOM MATRICES 



In the sequel, we will see how eigenvalues of some classical matrix- valued pro- 
cesses and radial Dunkl processes are interelated using SDE. This connection was 
already checked by physicists throughout eigenvalues probability densities and 
Fokker-Planck equations for parameter-dependent random matrices ([12]). As we 
mentioned in the introduction, the A^-i-type is connected to symmetric and Her- 
mitian Brownain motions. Set k := f3/2, f5 > 0, then such a process will be called 
/5-Dyson, referring to the Dyson model when (3 = 2. This parameter is called the 
Dyson index. Henceforth, we will adopt new notation for the eigenvalues process, 
we wil write A instead of X. 



5.1. The Bm-type.: /?-Laguerre processes. The Bm system turns out to be 
related to eigenvalues of Wishart and Laguerre processes which satisfy the following 
stochastic differential system (see [5], [13]): 



Xi{t) + Xk{t) 
^ \i{t) - Afc(t) 



dt 



1 < i < m. 



for (3 = 1,2 and 5 > m + l,m respectively, where (z/j)j are independent Brownian 
motions and Ai(0) > ■ ■ ■ > Am(0). Recall that the process remains strictly positive 
if it is initially strictly positive. This suggests to define the /3-Laguerre process as 
the solution, when it exists, of : 



d\i{t) = 2^/^t)dui{t)+p 



kf^i 



Ht) + Xkjt) 
kit) - Xk{t) 



dt 



1 < i < m. 



t < tAR 



0, 



where Rq = mi{t, Xm{t) = 0}, [3,6 > and with Ai(0) > ■ ■ • > Am(0) > 0. It is 
very easy to see that : 



dDt:=dll[Xi{t)] =2Dt 



,j=i 



for alH < r A -Rq where i? is a standard Brownian motion. It follows that Vr G 
diHD,)) = 2 E —-dB, + [(3{5 - m) + /? - 2] E jjwdt 

m _ 

(i(det(A)") = Mt + r[/3((5-m + l) + 2r-2]D[ V— -rft 

7^ 



where Mt = 2rDl^Y.T=i '^/\{t)dBt. From these two SDE, we can argue as m 
the Wishart and Laguerre cases that Rq > r a.s. when 6>m — 1 + 2/ j3 (choose 
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2r = 2 — j3{6 — m + 1) <0 when 6>m — 1 + 2/(3 then use McKean's argument). 
Set Tj := \/Xi, then, for t < r A Rq: 



dri{t) = dvi{t) + 
= dvi{t) + 



2r,(i) 



2ri{t) 



_ri(t)-r,(t) r,(t)+r,-(t) 



dvAt) H /^c/t + ki 



nit) -r jit) ri{t)+rj{t) 



dt 



with 2kQ = P{6 — m + 1) — 1, 2fci = /3. Consequently, the process r = (ri, . . . , r^) 
defined for all t < r A i?o is a i?m-radial Dunkl process. Using Theorem [H one 
claim that the SDE above has a unique strong solution for all t > and all (3, 6 
such that ko, ki > 0. This strengthen results from matrix theory : in the Wishart 
setting (P = 1), the strong uniqueness holds for 5 > m and in the Laguerre case 
{(3 = 2), it holds for 6 > m — 1/2. Besides, the generalized Bessel function is given 
by([2])0: 



77\ Yl Dk{x,wy) = 



[2 IP), [3 5 



\W , 

SO that ([T]) writes 

(9) p^'=^(x,y) = 



2 ' 2 ' 2 



-EE 

p=0 T 



2 )r 4'^'\Up\ 



\w\ 



Of: 2 2^ 



i=l 



where V stands for the Vandermonde function. Using the variable change y ^— y/y, 
the semi-group density of the /3-Laguerre process writes : 



koM . 



^7+2fci+m/2 



2 2t 2t 



l[{y^)'^o~l/2y2k,^y^dy 



i=l 



For X = and t = 1, we recover the same p.d.f. given in for /3-Laguerre 
ensemble. 

Remarks. 1/ Recall that for all a G -R, we set l{a) = k{a) — 1/2. Hence, in the 
i?m-case, Iq = ko — 1/2, /i = fci — 1/2. For — /, all corresponding parameters will be 
primed. For instance, —Iq = k'^ — 1/2, —li = k[ — 1/2. Let us consider a Wishart 
process of dimension S' such that m — l<5'<m + l ([5J), A;^ = 1/2 (/?' = 1) and 
k'^ = {5'-m)/2 => -h = 0, -lo = {6'-m-l)/2 < 0. Set 6' = m + l-2u, < z/ < 



With the same notations used in [2j, one has /?a'/2 = ko, a = ko — 1/2, (3 = 2/a ^ a + q 
P6/2. 
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1/2, then, = 0, /i = z/ ^ fci = 1/2 (/5 = 1) and fco = z/ + 1/2 (5 = m + 1 + 2v). 
Results of 4.1 writes: 



m / 2 \ ^ 



m + 1 5 a;^ , 
2 '2'2t- 



which fits the expression aheady derived in [15]. When k'Q = k[ = {—lo = —h = 
—1/2), then = ki = 1{(3 = 2, 5 = m + 1/2) and the Jack polynomials fits the 
Schur functions (see [311 )• ^^at case, the following representation holds ([21!) 



det(2;™ \^i{a -j + l,6-j + l, Zi)i<i^j<^ 



iFi'\a,b,z) 



Viz) 



where i^^i denotes the univariate hypergeometric function. Hence, the tail distri- 
bution writes : 



p-\To>t) = Ckdet 



2 \ m-j+1/2 



i iv 2 ^ ' 2 2t' 



The corresponding process is the Brownian motion in the Weyl chamber of 5-type. 
For Laguerre processes ([13j) of dimension 6 such that m— 1/2 < 6 < m, one should 
apply results derived in section 4.2. Take k[ = l{/3' = 2) and k'^ = 6' — m + 1/2 =^ 
li = 1/2, lo = 6' - m := -z/ with < z/ < 1/2. Thus ri = 1/2, = v ^ ki = 
1 (/9 = 2) and ko = v + 1/2 {5 = m + v) so that : 



m / 2 \ ^ 

i^i(To>t) = C.n(|) e-N^/-,F«(m,5, 

8=1 ^ ^ 



2t- 



2/Recall that when /3 = 2, o-^i^^ has a determinantal representation (see [Mj) 
yielding to Konig and O'Connell result on the l^-transform of m-independent 
squared Bessel processes (BESQs) constrained never to collide (or stay in the 
Am-i-type Weyl chamber, see [31]). Similar results holds for y4m_i-type root sys- 
tem with Brownian motions instead of BESQs. Nevertheless, when k^ = ki = 
1 {(3 = 2,S = m + 1/2), a similar interpretation involving m- independent Brown- 
ian motions killed when they reaches holds. However, the Vandermonde function 
may be replaced by the product over positive roots. In this case, the eigenvalues 
process is known as the BM in the Weyl chamber of type Bm (see [23] for further 
details and other root systems). Since 7 = and from ([21], [13]) : 



.Fi^)((m+a 



'i)l<i<ry 



{m+bj] 



l<i<S5 



det[r^s{{ai + l)i<i<r, (bj + l)i<j<s, Xiyf)]ij 



V{x)V{y) 
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transforms to : 



h(v) e-(l^l'+lyl')/2* ™ /X V- 
l,j=l 



det 



h{y) e-(l-f +l?/P)/2* 
— ^ 7—n\ det 



'h{x) 



fm/2 



3 {Xi'llj)'^^ "' 

2' 4^2 



where h is the product over positive roots. Besides, the following holds (see [6j) 



Thus, 



h{y) 1 
''h{x) (27rt)™/2 



sinh(2A/^). 



2v^ 



sinh 



det [iVt(?/j - Xi) - Ntivj + Xi)] 



where Nt{u) = (l/v^27rt)e "^/2*^ which fits Grabiner's result ([23] page 186). This 
is in agreement with the generator since Ah = ([23J) and 

m 

^f = Af + r(log h,f) = Af + Y, d,{log h)dj, 

1=1 

where T is the so-called "operateur du carre du champ" (see [2Z] Chap. VIII). 
Besides, for m = 1, r is a Bessel process of dimension 26 = 3 and the expression 
inside the determinant in the second line is exactly the semi-group of the Brownian 
motion killed when it reaches (see [37] p. 87). 



5.2. Generalized Bessel function in the Dm case. In the classification of root 
systems, the Am-i and Bm are known to be "irreducible" and both of them cor- 
respond to some matrix processes. Another one, yet with no underlying matrices, 
is the root system defined by (see ^| p. 42) 

R = {ie-t ± Cj, 1 < i < j < m}, i?+ = {cj ± ej, 1 < i < j < m} 

There is one conjugacy class so that k{a) = ki. Grabiner's result reads for the 
Brownian motion in the Weyl chamber of -Dm-type {ki = 1) : 

I, ^ _ ^(y^) dei[Nt{y^ - Xj) - Ntivi + xj)] + det[Nt{y^ - xj) + Nt{yi + x^)] 

Pt i-^) y) 



V{x^) 

Cm ^-(|^|2+|^|2y2^ det [smh(xiyj/t)] + det [cosh{x,yj/t)] ^^^^^2^ 



V{xyAt^)V{y^) 
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where 7 = m{m — 1). The second term in the sum involves the transition density 
of a reflected Brownian motion (|-B|, see [37] p. 81). A natural way to interpret 
the announced formula is that the Weyl chamber is given by : 

C = {X eR"", Xi> ■■■ > Xm^i > \Xm\}, 

SO that C fits the -B^-Weyl chamber when Xm > 0, otherwise, it is its con- 
juguate with respect to Se„ since this simple reflection acts only on Xm and retains 
the others. With the help of the determinantal formula used before ([21], [13]), 

0=^1(91 ^) 



Csinh(2^)/^ and o^i(l/2, z) = cosh(2A/i) ([0]), one writes: 



-(|x|2 + |y|2)/2t 



n 

.1=1 



~2t' 



(1) 



1 X2 

2' 2t' 2t ' ° ^ 



1 y"^ 

m , — , — 

2 ' 2t ' 2t 



With regard to ([T]) and setting g = 1 + (m — l)ki, it is natural to prove that 
Proposition 6. 



x'^ y^ 



I I w&W i=l \ / \ 

Proof, it relies on the following characterization ([M]): given a reduced root 



1 

2' 2t' 2t 



system R with finite reflection group W, D 



tion valued 1 at x = satisfying A^D^^x, ■) = \x\'^D 



is the unique l^-invariant func- 



w. 



It is easy to see 



that the function above is ly-invariant since W is the semi-direct product of the 
symmetric group Sm and (Z/2Z)'"~^ acting by even sign changes. However, it is 
not for the finite reflection group associated to the Bm root system due to the term 
multiplying the first hypergeometric series. In the Dm case, the Dunkl Laplacian 
writes on VT-invariant functions : 



i=l 



+ 

m - Vj Vi + Vj 



d,. 



Let : 



f{x,y) 

d{x,y) 



(lA-i) 



p(iA-i) 



q + 



2' 2 ' 2 
1 

2' Y' ~2 



l[{x^y^/2) 



i=l 



considered as functions of the variable y such that the generalized Bessel function 
is proportional to df + g and Afc[(i/ + g] = Ak{df) + Ak{g). Recall that: 



2 2 
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V\y' 



Note also that is a particular case of the Dunkl Laplacian considered for the 
Bm root system when ko = 0. As a result 



A,g{x,y) = Ai''-\ko = 0) 



1 



\W\ 



<x,x> g{x,y) 



For the remaining term, note that both d and / are H^-invariant. Write = 
Yl^i where Tj is the difference Dunkl operator (see [3H] P- 5). Then using (see 
[38j p. 6), one has the derivation formula Ti{df) = dTi{f) + fTi{d). It gives that 

m 

A,{df) = dAkif) + fA,{d) + 2 Y,{Tmmu)) 



i=l 



Moreover, Ti{d) = di{d) and Tj(/) = di{f) by ly-invariance. Next we compute : 

4:kid(x, y) \^ ^ = 



Ak{d){x, y) = 2kid{x, — 



1 1 
+ 

Vi - Vj Vi + Vj 



As a result : 

m 

A,(d/)(x,2/) = [ciA,(/)](x,2/) + 25^[(a,(rf))(a,(/))](x,i/) 



i=l 



m ^ 

[dA,{mx,y) + 2[dY,-id,{fm^,y) 

i=i y^ 

m ^ 

d{x,y)[Au + 2Y,-di]f{^.y) = d{x,y)Af-\k, = l)f{x,y) 



i=l 



When fco = 1, / fits the generalized Bessel function in the case =^ Ai:{df){x, y) 
{l/\W\) < x,x > df{x,y). Finally : 



Ak 



J2 Dk{x,wy) 



~ lie • ijC 



^ Dk{x,wy) 



6. Alcove-valued process 



6.1. /3-Jacobi processes. Recall that the eigenvalues of the real Jacobi matrix 
process of parameters (p, q) (see [iGj for facts on this process) satisfy : 



dXi{t) = 2V(Ai(t)(l - \i{t))dvi{t) + 



{p-{p + q)\i{t)) + Y, 



Xi{t){l-Xj{t)) + Xj{t){l-\iit)) 
Xi{t)-Xj{t) 



forO < Xm{0) < ■ ■ ■ < Ai(0) < 1 and alH < inf{s > 0, A^(s) = OorAi(s) = l}Ar. 
The /3- Jacobi process is defined as a solution, whenever it exists, of the SDK 
differring from the one above by a parameter /3 > in front of the bracket. It 
is easy to see that if A is a /3- Jacobi process of parameters (p, q), then 1 — A is a 
/3- Jacobi process of parameters {q,p)- As mentioned in the introductory part, the 
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connection with root systems is not new in its own ([3]) however we prefer giving 
some details of this transition. Setting Aj = sin^ </)j then = arcsinV%" := ■s(Aj) 
and < (j)m < ■ • ■ < 01 < 7r/2. The first and second derivatives of s are given by: 

1 , _ 2(2sin2 0,- 1) _ 2cos20i 

sm 2(pi sm 20j sm 20, 

Using 

sin^ 0j — sin^ 0j = 2 sin(0j + (pj) sin(0j — (pj) 



sin^ 0j cos^ 0j + cos^ (pi sin^ 0j = sin^(0j + (pj) + sin^(0j — (pj) 
then, Ito's formula gives: 

d<pAt) = dui{t) + ^(P:iip±iW^ _ cot 2^idt 

sin 20j 

P dt ^ sm^i(P,it) + (Pj{t)) + sm\(Pi{t)-(Pj{t)) 
2 sin20i(t) ^ sin(0i(t) + sin(0,(t) - 

Writing sin^*;/)^ = (1 — cos20^)/2, then 

^ du.it) + ;3 V^-^ + ^^^^^ cot 2MtyM 

(3 dt ^ ?,m\(P,{t) + (Pj{t)) + sm\(Pi{t)-(Pj{t)) 
2 sin20i(t) ^ sin(0,(t) + 0,(t)) sin(0i(t) - 

where < (pmifi) < ■ ■ ■ < 0i(O) < 7r/2. Moreover, 

sin20i = [cot(0j + (pj) + cot(0j — (pj)] sin(0j + (pj) sin(0j — (pj) 
which gives 

mt) - Mt) + ^^-^ + ^^^^ cot 2(p.m 

P^ [l/sm\(P,{t)+(P,m + [l/sin^(0,(t) 

2 ^ cot (0, (t) + 4>j{t)) + cot (0i (t) - 0,- (t) ) 
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Using 1 + cot^ z = \l sin^ z, then 

. + ^^-1^ + cot 2mii 

2 ^ cot (0, (t) + 0,- (t) ) + cot (0i (t ) - 0,- (t) ) 



E 

Using 

and 



2 sin20i(t) 

1 - cot(0,(t) + cot(0,(t) - ^ cot(0i(t) + + cot(0i(t) -0,-(t)) 

COt(0i(t) + + COt(0i(t) - 



+ 



dt 



— cot{u + v) 
1 2 cos^ (/),• 



1 — cot('u) cot(f ) 
cot('u) + cot(u) 



COS. 



cot — cot : 



i{t) = dVi{t) + 



sin 20j 2 sin 0j cos 

we finally obtain 
(10) 

/Cq cot 0j + ki cot 2(j)i{t)dt + /C2 ^^[cot(0j + 0j) + COt(0j — 

where 

(11) 2ko = Pip-q), fci = /3(g-(m-l))-l, 2^2 = /5. 

Easy computations show that 7r/2 — satisfies ffTUj) with (p, g) intertwined. 

6.2. Eigenfunctions and Heckman-Opdam's functions. Let k2 > and =Sf 

be the generator of 0, then the eigenfunctions of ^ are given by Gauss hyperge- 
ometric series: in fact, let ^ be the generator of (Ai, . . . , A^) (see [16] p. 135 for 

(3 = 1): 



dt 



1=1 



i=l 



P 



Xi — Xi 



2 J2 ^^(1 - ^^)^■ + - (m - 1) - (p + g - 2(m - 1))A,]9, + 2/? ^ 



1=1 



1=1 



a; 



Ai(l - A, 

A,; — Ao 



■ft 



From Equation (IHl) (A;2 plays the role of fci), one can see that 2Fi^^''^\a,b,c; A) is 
the unique symmetric analytic function u such that u{0) = 1 which satisfies 

£/u{X) = 2mahu{X), 2c = (3p = 2ko + ki + 2k2{m-l) + 1, 2a + 26+1- 2c = k^. 
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with ki, < i < 2 cited in ffTTj) . Setting sin^ := (sin^ 01, ... , sin^ 0^), then j^/ 
transforms to Hence: 

=^[M(sin^ 0)] = 2mab [^(sin^ 0)] 

In the same spirit, one can also interpret ^ as the "radial part" of the trigonomet- 
ric version Dunkl-Cherednik Laplacian (with cot replacing coth, [3], [36]). By radial 
part, we mean the restriction on 14^- invariant functions. Besides, this Laplacian 
arises, as for Dunkl and Cherednik-Dunkl ones, from differential-difference first- 
order operators. However, this comes beyond the spirit of this work and will not 
be done here. 

6.3. Existence and uniqueness of a strong solution. The involved root sys- 
tem is the non reduced BCm defined by 

R = {icj, ±2ei, 1 < i < m, ±(ej ± ej), 1 < i < j < m} 

-R+ = {cj, 2ei, 1 < i < m, (cj ± Cj), I < i < j < m} 

A = {ci - Ci+i, 1 < i < m - 1, Cm} 

When ko = 0{p = q), it reduces to the reduced Cm system 

R = {icj ± Cj, 1 < 2 < j < m, ±2ej, 1 < i < m} 

R+ = {cj ± Cj, 1 < z < j < m, 2ej, 1 < i < m} 

A = {cj - Cj+i, 1 < i < m - 1, 2em} 

and it is known as the ultraspheric case. The Weyl group action on gives rise 
to three orbits so that the multiplicity function is given hj k = (fco, ki, ^2). Setting 
0j := 0i/7r, then the process is valued in the positive Weyl alcove (see pTj) defined 
by: 

i = {0 G M", < a,0 >> OVa G A < a,0 >< 1} 

where a = 2ei is the highest positive root (that is a — a G i?+ Va G R, see \27\). 
The associated affine Weyl group Wa is the semi-direct product of W and the 
translation group corresponding to the coroot lattice (Z-span of {2q:/||q:|P, a G 
R})- The generator writes in this case: 

^^(0):=^A^?(0)-< V^7(0),V<I>(0)>, $(0) = - A:(«)logsin(<«,0>) 

Thus, with minor modifications. Theorem [T] states that ffTUj) has a unique strong 
solution for all t > subject to ko > 0, ki > 0, k2 > P > 0, p > q > 
{m-l) + l/p. Applying this to vr/2-0, this holds ioi /3 > 0, q > p > {■m-l) + l//3. 
Since the ultraspheric case still involves a root system, then ( ITOl) has a unique 
strong solution for p Aq > {m — 1) + 1/(3 which simplifies to p A g > m in the real 
case /? = 1 and pAq>m — l/2m the complex one (3 = 2. Theorem [1] is modified 
as follows: dA = U^gA-f^a U H^^i where 

Ha,l = {05 < «, >= 1} = {0, TT— < (5, >= 0} 
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Compared with (jH), the convex function x ^ — ln(< a, x >) should be substituted 
by (/) I— — ln(sin(< a, >)) and one has to deal with an additional term in the 
expression of the boundary process {Lt)t>o '■ l{7r-<«,(^>=o}- Then the occupation 
density formula writes: 

/ U^ij- < a,(t) >)\6\a)\da =< a,a > / \6'{ti- < a, Xs >)\ds 
Jo Jo 

[ \e'{<a,X,>)\ds 
Jo 



--< a,a > 



since cot(7r — z) = — cot(z). Hence, the same proof applies and Lemma 1 remains 
valid for a G -R+. Besides, either it will exist a G A such that < a,x >= and 
Lemma [2] applies, or we will need to prove that < n{x), a >^ if x belongs only 
to Ha,i- Let us first recall that the highest root is the unique positive root such 
that a — a G R+ for all a G -R+. Thus it may be written as a = J2aeA ^a'^ where 
tta > 1- Else, if there exists ao E A such that < 1 and since a must be greater 
than all simple roots (in particular greater than oq) then 

a — ao = {ttaQ — l)tto + o,aa = Cagao + Caa 

for some Cq, > 0. Our claim follows from the fact that A is a basis. Next, it is not 
difficult to see from the definition of n{x) and the fact that < a,x » for all 
a G A that n{x) is colinear to — XloeA follows that 

< n{x), a >= — c < a,a >= —c aa < a,9 > 

qGA QGA6»eA 

If < n{x), a >= 0, then 

||^a||^ = ^^<a,6' >< ^^Oq, < a,6' >= 

aeA ogA 0gA cjGA eGA 

which implies that n{x) = 0. ■ 

6.4. Brownian motion in Weyl alcoves. Let 

hi{(j)) := JJ^ sin(<a,0>) 

OG-R+ 

Then, hi is strictly positive on A and vanishes for cj) G dA. One can also show 
that (l/2)A/ii = ch for some strictly negative constant c. Let P^^ denote the semi 
group given by : 
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where Pt denotes the semi group of the process consisting of m-independent BMs 
in A killed when it first reaches dA. The corresponding generator writes : 



in' . m 



i=l 



which fits our generator for k2 = 1 {l3 = 2) , ki = 2 [q = m + 1/2), fco = l(p = 
q + 1 = m + 3/2). In the ultraspheric case, this becomes P = 2, p = q = m + 1/2. 
In both cases, these parameters correspond to the process consisting of m BMs 
constrained to stay in the BCm and Cm- Weyl alcoves respectively. Note that 
p, q are not integers which means that these processes BM can not be realized as 
eigenvalues processes of complex matrix Jacobi processes which is also the case for 
the BM in the -B^-Weyl chamber since 6 = m + 1/2. 

6.5. The first hitting time Tq. We define similarly the first hitting time of 
alcove's walls by Tq = mf{t > 0, {(j)(t)/n) G dA} = A mf{Ta, a e A}, where 

fa := inf {t > 0, < a, 0(t) >= 0}, 

fa := inf {t > 0, < a, 0(t) >= 20i = vr}, 

and (j) is the unique strong solution for alH > of 0: 

d(j){t) = du{t) + k{a) cot(< a, 0(t) >)a dt, G A, 

for the non reduced root system R = BCm with k{a) > for all a and p A q > 
{m — 1) + 1/p. Let us focus on Tq^ for some ao G A. We shall distinguish two 
cases : 

6.5.1. ao = Ci — Cj+i, 1 < i < m — 1. The same scheme described in the proof 
of Proposition [T] applies here since the main ingredients used there are the SDE 
and the fact that cTo(a) G R+ \i a ^ a^. The second assertion follows from 
cro(2ej) = 2o-o(ej) = 2(5ijei+i + 5(i+i)jei + l{jyijyi+i}ej) G -R+. As a result, one 
writes for alH > 0: 

d < ao, 4>(t) >= I |ao| Mt* + ^2! |«o| 1^ cot < ao, 4>(t) > dt + k{a)a{a) cot < a, (j)(t) > dt 

where a{a) =< ao,a >. 

d < ao, (pit) >= \ \ao\ \d'jt + ^2! Ictol P cot(< ao, 4>{t) >)dt + Ft 



where 



Ft = /i;(a)a(Q;)[cot(< a, 0(t) >) — cot(< (To(a), 0(t) >)], 



aeR+\ao 
a(Q)>0 



^k{2e^) = fci/2 for all 1 < i < m. 
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where (Tq = (Tao- This drift is strictly negative on {T^g = 00} since (f) 1-^ cot cf) is a 
decreasing function, < ao,(j)(t) » and since : 

< a - aoia), (pit) >= 2-^^^ < ao, (pit) » 0. 

I I'^l 1 

This imphes that W^(it > 0, < ao, (f){t) >< Zt) = 1 where 0(0) = x and : 

dZt = \ \ao\\d'yt + | |q;o| ^^2 cotiZt)dt, Zq =< ao, 0(0) >= x 

on the same probabihty space. Using ( |TOl) with (3 = l,m = 1, one can easily see 
that {Zt)t>o = (arcsin t>o where J is a one dimensional Jacobi process of 

parameters d = 2k2 + l,d' = 1 (see ^|) : that is : 

dJt = 2^Jt{l - Jt)d-ft + {d-{d+ l)Jt)dt, < fca < 1/2 ^ < < 2. 

As J hits almost surely when < < 2 (use the skew product in [43j and 
properties of squared Bessel processes), then so does Z and by the way < ao, (p > 
for /c2 < 1/2 ^ < 00 a. s. 

6.5.2. ao = Cm- Compared with the previous case, the difference arises from the 
fact that cro(a) G R+ if a E R+\{em, 2em} and the latter is easily checked since for 
a = Cj ± ej this amounts to consider the reduced root system Bm, else for a = Cj 
with i ^ m, (To(ej) = Cj. According to this, one gets : 

d < ao, (pit) >= d(pmit) = d'jt + ko cot(0m(t))(it + ki cot(20m) + Ft 

where 

Ft= ^2 fc(a)a(a) [cot(< a, (pit) >) - cot(< croia),(pit) >)] 

aeB.+ \{em ,2em} 
a{a)>0 

where R\_ = {ei — Cj, 1 < i < j < m}. Using once again (fTUj) . we shall compare 
this process with (arcsin \/Jt)t>o where 

dJt = 2^ J til- Jt)d-it + id-id + d')Jt)dt, d' = ki + l,d = 2ko + ki + l. 

Hence, f^^ < 00 a.s. if < 2A;o + /ci < 1/2 /3p - (/5(m - 1) < 2. This agrees 
with the case m = 1 for which p < 2 (use the skew product in |13]). Finally, note 
that since a(a) = for a G {e,, 2ej, 1 < i < m — 1}, F only involves k2 = (3 which 
is independent from p, q. Keeping in mind that 7r/2 — is still a /5- Jacobi process 
with ip, q) intertwined which has no effect on the strict negativity of F by the 
above remark, we conclude that Tq, < cxd for < /3g — /3(m — 1) < 2. ■ 
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6.6. Semi-group density. We end this paper by giving the semi group density 
of the /3-Jacobi process. Before proceeding, we briefly consider two cases for which 
we can write down the semi-group density: the univariate case and the complex 
Hermitian one {(3 = 2). Let P^'* denote the Jacobi polynomial of degree n defined 

by my- 

pr,s^X) := ^^^^2^1 (-n, n + r + . + 1, r + 1; , 

for A G [— 1, 1], r, s > — 1, where 2-^1 is the univariate Gauss hypergeometric func- 
tion. These polynomials are orthogonal with respect to the measure Z^''^{X)dX : = 
(1 — A)''(l + Xj'^dX and the associated inner product in L^([— 1, 1]) given by 

</,^?>i.([_i,i]):= / f{X)g{X)Z''^{X)dX 

Moreover, (P^''*)„>o form a complete set of this Hilbert space and satisfy 

J^Vl^dl + [{s-r)-{s + r + 2)A]9a} P^ (A) = -n{n + r + s + 1)P:'^(A) 

The above eigenoperator defines a diffusion which is related to the one we con- 
sidered with m = 1 via the map A 1-^ (1 — A)/2 and a deterministic time change 
(t I— > t/2). The semi group density w.r.t Lebesgue measure is written (see 

00 

pt'io,x) = Y,e-'^--'p::'%e)p:^'%x)w'''%x) 

n=0 

where r„ denotes the eigenvalues above, (-P^''')n are orthonormal polynomials, p = 
2(r + l),q = 2(s -|- 1) and W^''^{X)dX is the probability measure corresponding 
to the measure Z^''^{X)dX. No closed forms seems to be known for this density, 
nonetheless an attempt to get a handier expression was tried in [13] . Multivariate 
analogs appeared in literature ([S], [2H], [32] for instance) and are obtained by 
applying the Gram-Schmidt orthogonalization to the symmetric Jack polynomials 
w.r.t. measure 

m 

Z:;:'^{X)dX:=l[Xl{l-Xir H \K-X/dX,...dXm 

i=l l<«<i<»7i 

We shall denote them0 by P^'**'^ for a given partition r (instead of G"'^ used in 
literature) and stress that some of the properties cited above extend to the higher 
dimensional case ([32]): an expansion in terms of 2-^i^^^''(— b, c. A) exists for r = 
(Z'") with m components all equal to /; (P^'*'^), where r is a partition of length 
< m, form a basis of the Hilbert space ^^([0, l]'", W;;i'^^{X)dX) where W^''^{X)dX 
is the normalization Z^'^{X)dX in order to be a probability measure ([32])- The 
normalizing constant is given by a McDonald- Selberg integral computed in 



The normalization is different from the one used in both 3J and [32j . 
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Moreover, (P^'*'^),- are the unique symmetric polynomial eigenfunctions of the 
Laplace Beltrami operator (thus defined on [0, 1]™) with — {m — 1)) = 
2(r + 1), — (m — 1)) = 2(s + 1), associated with the eigenvalues 



(12) 2r; 



5^ r,(r, - 1 - Pit - 1)) + \r\{r + s + (3{m - 1) + 2) 



r = n. 



However, with regard to the strong uniqueness for all t > previously derived, we 
shall restrict ourselves to p /\q> (m — 1) + [3{q — (m — 1)) > 1 is equivalent 
to s > —1/2 and f3{p — (m — 1)) > 1 is equivalent to r > —1/2. As a result, 
r,s > -1/2. 

It is known that the eigenvalues process of the complex Hermitian Jacobi process 
(or 2- Jacobi process) is the /i-transform (in the Doob sense) ioi h = V ol a, process 
whose components are real Jacobi processes of parameters 2(p — (m — 1)) = 2(r + 
1), 2(g — (m — 1)) = 2(s + 1) constrained to never collide (or to stay in the A^-i- 
type Weyl chamber). Here, V denotes as usual the Vandermonde function. More 
precisely, V is an eigenfunction of the generator of the one dimensional Jacobi 
process of parameters (p, g) (see appendix in [16]), say L, that is 



I.V = cV 



-m{m — 1) 



2(m — 2) p + q 



V 



Noting that the parameters r, s are the same both in the univariate and in the 
multivariate cases, it follows by Karlin-McGregor formula ([2S]) that the semi 
group density writes on{0<Am<---<Ai<l} 



r,s,2/ 



A) :=e- 



V{0) 



det K]]e-2"("+^+^+i)*P;'^(^i)P;'^(Aj)Vr^'^(Aj 



vn=0 



.n=Q 



e-^'det I 5]]e-2"("+^+^+i)*P„^'^(^,)P„^'^(Aj 

E 



v{e)v{X) 



.crieS'm ni,...,nm>0 



-^'ETLiMrii+r+s- 



o-i(j)y 



i=l 



-ct 



.(Tl,(T2G5m ni>--->nm>0 



E 



-2EII1 »i<T2(«) ('^<T2(i)+»'+S+l) 



v{e)v{X) 



Note that, for a given partition (ni > ■ ■ • > rim > 0) and a permutation (J2 G S'„ 
one has 

m m 

5]^'^<72(i)(^<T2(i) + r + s + l) = ^nj(ni + r + s + 1) 

i=l i=l 
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W^^'^(A) 
^W^(A) 



Thus summing first over (Xi with the change of variables a = (TiO'2, one gets: 

' ^ ' ^ ...k.o " no) vix) 

,det[P„Y(^,)],,detKf(A,)],, 



ni>-->nm>0 

Set rij = Tj + m — z, then ri >> ■ ■ ■ > > 0. Moreover, with regard to ffT^ . one 
easily check that 



^2 



ri(ri + r + s + 1 + 2(m - i)) 



i=l 

SO that 

m 

^nj(ni + r + s + 1) = r^_^ - c/2 

i=l 

The final result writes 



;<-''.S'2//i ^-2ra ^ det[P^;^^_^(6'j)]j.j det[P^;^^_^(Aj)]ij 2^ ^ ^ 
K, {e,X)= e (A) 

Tl>...rm>0 

where we used the determinantal representation of the Jacobi multivariate poly- 
nomials in the complex cas^ (see [32]) : 

pr.,s,2.^^ _ det[P;;+^^,(A,)],,, 

From these observations, it is natural to claim that : 

Proposition 7. The semi group density of the (3- Jacobi process is given by 

00 

(13) Kl'-^^e, X):=Y.J2 e-^"-*^;'^'^(^)^;'^'"(A)W^r(A)l{o<A„<...<A,<i} 

n=0 |r|=n 

with respect to dX. As a result, it is positive. 

Proof, given a bounded symmetric function / on [0, l]*", define 
„ 00 
TJie) := / /(A) E e-''^"-*P;'^'^(^)P;'^'^(A)W^^^(A)rfA 

^0<A„<-<Ai<l „=0|r|=n 

for = (0 < < ■ ■ • < < 1) and Tq/ = /. The above expression makes sense: 
this uses the boundness of /, the exponential term with strictly positive t and 
Fubini Theorem. Besides, T^l = 1 and ||Ti|| is bounded for all t > 0. The first 



We adopt a different normalization since we consider orthonormal polynomials. 
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claim follows easily from the orthogonality of the Jacobi polynomials and Pq = 1 
so that the only non zero term is that correponding to n = 0. The second one 
is obvious for t = and uses the exponential term when t > e > 0. One also 
easily checks that TfT^ = Tt+s and that ^Ttf{\) = dtTtf\X) using the dominated 
convergence theorem. Now, let us consider the Cauchy problem associated to =Sf : 

^(t,A)=ifn,(t,A) 

«/(o,-) = /, 

where m/ G Ci'2(M* x 5 := {0 < A„ < ■ ■ • < Ai < 1}) na(M+ nS) with reflecting 
boundary condition : 

< Vn(t,A),r2(A) >= 0, (t, A) G M+ x 95 

where n(A) is a unitary inward normal vector at Define Ut{f){X) := Uf(t,X). It 
is shown ([12]) that there is a unique solution to the Cauchy problem with ini- 
tial condition. Consequently, {Tt)t>o is the semi group of the eigenvalues process 
{X{t))t>o with density given by Kl'^'^ . ■ 

Remark. As the reader can check, the computations performed in the complex 
Hermitian case do not restrict to Jacobi polynomials. We only used the determi- 
nantal representation in terms of their univariate counterparts. As a result, one 
gets similar formulas replacing Jacobi by Hermite and Laguerre polynomials. 

Now, we are able to answer some open questions left in (TSj . For the real Jacobi 
matrix (/? = 1), it is known that for p /\ q > m — 1 and if the eigenvalues are 
distinct at time t = 0, then they remain distinct forever. It is then natural to 
wonder if this remains valid when starting from non distinct eigenvalues (see [l6] 
p. 138-139). The Markov property together with the previous result for distinct 
eigenvalues are sufficient to claim that this is true provided that the eigenvalues 
semi group has a density which is absolutely continuous with respect to Lebesgue 
measure on R™. By virtue of Kl'^'^{6, 0), for p A q > m, 

PA(o)(Vt > 0, Vz ^ J, A,(t) ^ A,(t)) = 1, Ai(0) > • • ■ > A„(0). 

We argue in the same way to claim that for p A q > m + 1, the process will never 
hit the boundaries (0 and 1 for A or and tt/2 for 0) even if it did at time t = 0. 
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